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(57) Abstract 

A system and method for identifying characteristics of a physical system (10) applies substantially identical drive ensembles (I 15. 
1 17) to the physical system (10) and obtains corresponding responses (21) from the physical system (10). A repeatable bandwidth of the 
physical system (10) is estimated as a function of the applied drive ensembles (1 15, 1 17) and the corresponding obtained responses (21). 
Instructions can be provided on a computer readable medium to perform the method. 
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MEASURING SYSTEM REPEATABLE BANDWIDTH 
FOR SIMULATION TESTING 

CROSS-REFERENCE TO RELATED APPLICATION 
This application is based on and claims priority 
5 of U.S. provisional patent application serial number 
60/108,295 filed November 13, 1998, the content of 
which is hereby incorporated by reference in its 
entirety. Co-pending and commonly assigned U.S. patent 
application serial number 09/234,998 filed January 21, 
10 1999 is also incorporated by reference in their 
entirety . 

BACKGROUND OF THE INVENTION 
The present invention relates to control of 
•a system, machine or process, that is repetitive in 
15 nature or is amenable to at least some degree of 
rehearsal. More particularly, the present invention 
relates to determining a repeatable control bandwidth 
of. a vibration system, to or within a metric 
appropriate for the application. 
20 . Vibration systems that are capable of 

simulating loads and/or. motions applied to test 
specimens are generally known. Vibration systems are 
widely used for performance evaluation, durability 
tests, and various other purposes as they are highly 
25 effective in the development of products. For 
instance, it is quite common in the development of 
automobiles, motorcycles, or the like, to subject the 
vehicle or a substructure thereof to a laboratory 
environment that simulates operating environments such 
30 as a road or test track. Physical simulation in the 
laboratory involves a well-known method of data 
acquisition and analysis in order to develop drive 
signals that can be applied to the vibration system to 
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reproduce the operating environment . This method 
includes instrumenting the vehicle with transducers 
"remote" to the physical inputs of the., operating 
environment. Common remote transducers include, but 
5 are not limited to, strain gauges, accelerometers, and 
displacement sensors, which implicitly define the 
operating environment of interest. The vehicle is then 
driven in the same operating environment, while remote 
transducer responses., (internal loads and/or motions) 
0 are recorded to represent the "desired" response for 
the simulation. During simulation with the vehicle 
mounted to the vibration system, actuators of the 
vibration system are driven so as to reproduce the 
recorded remote transducer responses on the vehicle in 
5 the laboratory thereby replicating the desired 
response. 

However, before simulated testing can occur, 
the relationship between the input drive signals to 
the vibration system and the responses of the remote 
) transducers must be characterized in the laboratory. 
Typically, this "system identification" procedure 
involves obtaining a respective system model or 
transfer f unct ion of the complete physical system 
(e.g. vibration system, test specimen, and remote 
5 transducers) hereinafter referred to as the "physical 
system" . The . inverse of the system model is used to 
iteratively obtain suitable drive signals for the 
vibration system to achieve substantially the same, 
response from the remote transducers on the test 
specimen in ..the laboratory situation as was found in 
the operating environment . The iterative process can 
involve, for example, various methods of adjusting the 
drive signals iteratively until the response achieved 
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f rom the physical system is acceptably close to the 
desired response. 

As those skilled in the art would 
appreciate, this process of obtaining suitable drive 
5 signals is not altered when the remote transducers are 
not physically remote from the test system inputs 
(e.g. the case where "remote'' transducers are the 
feedback variables, such as . force or. motion, of the 
vibration system controller) . 
10 - Although the above-described system and 

method for obtaining drive signals for a vibration 
system has enjoyed substantial success, there is a 
continuing need to improve such systems. For example, 
a fundamental, limitation on the accuracy with which 
15 the desired operating responses can' be reproduced in 
the simulation test is the repeatability of the 
response of the physical system, as measured by the 
remote transducers, . to the same input drive signal 
(repeated) . Frequently, physical systems are only 
2 0 repeatably controllable to Within * an appropriate 
metric of accuracy over a limited frequency range. 
This limited frequency range, referred to herein as a 
"repeatable bandwidth", is the frequency range over 
which the system can be controlled with some * desired 
25 or necessary measure of repeatability. While the 
physical system will generally be controllable over 
larger frequency ranges, physical system 

characteristics result in repeatability exceeding the 
desired metric outside of the repeatable bandwidth. 
30 As discussed, one of the primary difficulties in 
performing laboratory simulation is to determine the 
frequency range over which the . simulation may_ be 
expected to be repeatably accurate. A common method of 
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predicting a simulation range is to use ordinary, 
Partial and multiple coherences whose results often do 
not correlate with the results obtained during the 
iterative process. Moreover, coherence measurements 
frequently, do not -provide- any recourse for improving 
the simulation bandwidth. 

Commonly, when developing drive signals for 
the physical system, an assumption is made that, the 
system will be used, and is repeatably controllable, 
0 over a particular frequency range. Considerable effort 
may be expended trying to achieve a drive signal that 
accurately reproduces the desired response over a 
frequency range that ' in fact may not be possible. 
Consequently, there is a. need to accurately estimate 
the repeatable bandwidth of the physical system prior 
to modeling and/or iteratively ' obtaining drive signals 
during the system identification phase.. 

SUMMARY OF THE INVENTION 
A system and method for' identifying 
characteristics of a physical system applies 
substantially identical drive ensembles to the 
physical system and obtains corresponding responses 
from the physical system. A repeatable bandwidth of 
the physical system is estimated as a function of the 
applied drive ensembles and the corresponding obtained 
responses. Instructions can be provided on a computer 
readable medium to perform the method. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of an exemplary 
environment for practicing the present invention. 

FIG. 2 is a computer for implementing the 
present invention. 

.FIG. 3 is a flpw chart illustrating steps 
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involved in an identification phase of a method of 
vibration testing. 

FIG. 4 is a diagrammatic illustration of a 
prior art technique for generating frames of a test 
5 drive signal. 

FIGS. 5 and 6 are diagrammatic illustrations 
of a first technique for generating frames of a test 
drive signal in accordance some embodiments of the 
. present invention. 
10 FIG. 7 is . a diagrammatic illustration of a 

second technique for generating frames of a test drive 
signal in accordance some embodiments, of the present 
invention. 

FIGS. 8 and 9 are plots illustrating an 
.15 aspect of calculating physical system repeatable 
bandwidth in accordance with the present invention. 
. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

FIG. 1 illustrates a physical system 10 . . The 
physical system 10 generally includes a vibration 
20 system 13 comprising a servo controller 14 and an 
actuator 15. In the schematic illustration of FIG. 1, 
the actuator 15 represents one or more actuators that 
are coupled through a suitable mechanical interface 16 
to a test specimen 18 . The servo controller 14 
25 provides an actuator command signal 19 to. the actuator 
15, which in turn, excites the test specimen 18. 
Suitable feedback 15A is provided from the actuator 15 
. to the servo controller 14 . One or more remote 
transducers 20 on the test specimen 18, such as 
30. displacement sensors, strain gauges, accelerometers , 
or the like, provide a measured. or actual response. 21. 

A physical system controller 23 receives the 
actual response 21 as feedback to compute a drive 17 
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as input to the physical system 10. In an embodiment 
in which the drive 17 is generated using an iterative 
process, the physical system controller 23 generates 
the drive 17 for the physical system 10 based on the 
5 .comparison of a desired response provided at 22 and 
the actual response 21 of the remote transducer 20 on 
the test specimen 18. Although illustrated in FIG. 1 
for. the single channel case, multiple channel 
embodiments with response 21 comprising N response 
10 components and the drive 17 comprising M drive 
components are typical and considered another, 
embodiment of the present invention. 

Although described herein where the physical 
system comprises the vibration system 13 and. remote 
15 transducer 20, aspects of the present invention 
described below can be applied to other physical 
systems. For instance, in a manufacturing process, the 
physical system includes the manufacturing machines 
(e.g. presses, molding apparatus, forming machines, 
etc.) and the drive 17 .provides command signals to 
said machines, and the actual response 21 comprises 
manual or automatic measured parameters of the 
manufactured article such as a critical dimension. 

FIG. 2 and the related discussion provide a 
25 brief, general description of a suitable computing 
environment in which the invention may be implemented. 
Although not required, the physical system controller 
23 will be described, at least in part, in the general 
context of computer-executable instructions, such as 
30 program modules, being executed by a computer 30. 
Generally, program modules include routine programs, 
objects, components,, data structures, etc., which 
perform particular tasks . or implement particular 



20 
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abstract data. types. The program modules are 
illustrated below using block diagrams and flowcharts. 
Those skilled in the art can implement the block 
diagrams and flowcharts to computer- executable 
5 instructions. Moreover, those skilled in the art will 
appreciate that the invention may be practiced with 
other computer system configurations, including .multi- 
processor systems, networked personal computers, mini 
computers/ main frame computers, and the like. The 
10 invention may also be practiced in distributed 
computing environments where tasks- are performed by 
remote processing devices that are linked through a 
communications network. In a distributed computer 
environment, program modules may be located in^ both 
15 local, and remote memory storage devices . 

The computer 3 0 illustrated in FIG. 2 
comprises a conventional personal or desktop computer 
having a central processing unit (CPU) 32, memory 34 
and a system bus. 36, which couples various system 
20 components, including the memory 34 to the CPU 32. The 
system bus 36 "may be any of several types of bus 
structures including a memory bus or a memory 
controller, a peripheral bus, and a local bus using 
any of a variety of bus architectures. The memory 34 
2 5 includes read only memory (ROM) and random access 
memory (RAM) . A basic input/output (BIOS) containing 
the basic routine that helps to transfer information 
between elements' within the computer 30, such as 
during start-up, is stored in ROM. Storage devices 38, 
30 such as a hard disk, a floppy disk drive, an optical 
disk drive, etc., are coupled to the system bus 3 6 and 
are used for storage of programs and data. It should 
be appreciated by those skilled in the art that other 
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types of computer readable media that are accessible 
by a computer, such' as magnetic cassettes, flash 
memory cards, digital video disks, random access 
memories, read only memories, and the like, may also 
5 be used as storage devices. Commonly, programs are 
loaded into memory 34 from at least one of the .storage 
devices 3 8 with or without accompanying data. 

An input device . 4 0 such as a keyboard, 
pointing device (mouse) , or the like, allows the user 
10 to provide commands to the computer. 30. A monitor 42 
or other type of output device is further connected to 
the system bus 36 via a suitable interface and 
provides feedback to the user. The desired response 22 
can be provided as an . input , to the computer 3 0 through 
15 a communications link, such as a modem, or through the 
removable media of the storage devices 38. The drive 
signals 17 . are provided to the physical system 10. of 
. FIG. 1 based on program modules executed by the 
computer 30 . and through a suitable . interface 44 
20 coupling the computer 30 to the vibration system 13. 
The interface 44 also receives the actual response 21. 
Using computer executed program modules, discussed 
below as methods, computer 3 0 calculates a repeatable 
bandwidth or frequency range of the physical system. 
25 The calculated repeatable bandwidth is then used in 
the modeling and/or drive signal generation processes. 

Referring now to FIG . 3 , the present 
invention can be described in conjunction with known 
methods for modeling the physical system 10 and 
3 0 obtaining the drive 17 to be applied thereto. Although 
described below with respect to a test vehicle, it 
should be understood that these known methods and the 
present invention . as illustrated in FIG. 3 are not 
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confined to testing only vehicles, but can be used on 
other types of test specimens and substructures or 
components thereof. In addition the system model can 
take many ■ forms including but not limited to linear 
5 and non-linear represented in frequency, time, state 
space or neural networks. Furthermore, the description 
is done assuming spectral analysis, based modeling 
estimation and implementation though operations can be 
carried by several other mathematical techniques (e.g. 
10 Adaptive Inverse Control (AIC) type models, parametric 
regression techniques such as Auto Regressive 
Exogenous (ARX) and State Space types of models, or 
combinations thereof) . 

Referring to FIG. . 3- at step 52, the. test 
15 vehicle is instrumented with the remote transducers 
20. At step 54, the vehicle is subjected to the field 
operating environment of interest and the .remote 
transducer responses are measured and recorded. For 
instance > the vehicle can be driven on a road or test 
2 0 track. The measured remote transducer responses, 
typically analog, are stored in the computer 3 0 in a 
digital format through analog- to-digital converters, 
as. is commonly known. 

Next, in a system identification phase 55, 
2 5 the bandwidth and the input /output model of the 
physical system 10 are determined. At step 56, 
modeling drive ensembles are applied to physical 
system 10 to obtain responses from the physical 
system At this point, it is helpful to provide 
30 definitions for terminology used, to describe the 
present invention. It must be noted, however, that 
slightly different definitions for the following 
terminology can be used without departing from the 
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scope of the present invention. For purposes of 
generating a model of the physical system, it is 
common to provide, on each channel, a "random" noise 
signal of a known duration sometimes referred to as a 
5 "frame". An exemplary duration for a frame is five 
seconds. The "random" noise signal is typically a 
combination of multiple periodic waveforms. An 
"ensemble" can be defined as a collection of 
individual channel signals over a frame or frames,. 
10 i.e. a subset of a time history and potentially the 
whole time history. For purposes of describing the 
invention, the phrase "modeling drive ensemble" can be 
defined as an ensemble used for system 'identification . 
However, a definition of the term "ensemble" being 
15 • equivalent to one or more frames of data for a single 
channel time history is also consistent with the 
present invention. 

At step 56 illustrated in FIG. 3, drive 
ensembles are applied to physical system 10 to obtain 
.20 a response. A method of applying th£ drive ensembles 
to the physical system in accordance with the first 
aspect of the present invention is discussed later 
with reference to FIGS. 5-7. With the exception of the 
method of applying the 'drive ensembles to the physical 
2 5 system in accordance with the first aspect of the 
present invention, step 56 of applying drive ensembles 
to the physical system to obtain the response of the 
physical system is similar to conventional methods of 
formulating a model of the physical system. 
30 Using the drive ensembles and the responses 

obtained from the physical system in response to 
application of the drive ensembles, an estimate of the 
repeatable bandwidth of the physical system is 
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calculated as illustrated generally at step 57. In 
accordance with the second aspect of the. present 
invention, system controller 23 (FIG. 1) calculates 
the . system repeatable bandwidth as a function of the 
difference between responses of physical system 10 to 
multiple replications of each drive ensemble. Due to 
characteristics of the physical system, different 
responses will be received during repeated 
applications of identical drive ensembles. as 
appreciated by those skilled in the art. Adding 
additional replications of each drive ensemble, in 
accordance with the first aspect of the . present 
invention, facilitates analysis of the system response 
to these ensembles to. calculate an estimate the system. 
15 repeatable bandwidth. This aspect of the present 
invention is described below in greater detail with 1 
reference to FIGS. 5-7. 

At step 58, a system model of physical 
system 10 is calculated. The physical system model can 
2 0 be calculated using known procedures. These procedures 
can utilize the responses to the drive ensembles 
applied to the physical system 10 at step 56. In other 
embodiments, the. system model and the system 
. repeatable bandwidth may be derived from a separate 
25 and potentially different set of drive ensembles. 

At step 58, an estimate of the model of the 
physical system 10 is calculated based on the input 
-drive applied and the remote transducer response 
obtained. In one embodiment, this is commonly known as 
3 0 the "frequency response function" (FRF) . 

Mathematically, the FRF is a N x M matrix wherein each 
element is a frequency dependent complex variable 
(gain and phase versus frequency) . The columns of the 
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matrix correspond to the inputs, while the rows 
correspond to the outputs . As appreciated by those 
skilled in the art, the FRF may also be obtained 
directly from prior tests using the physical system 10 
5 or other systems substantially similar to the physical 
system 10 - 

An inverse of the system model H(f)" 1 is 
needed to determine the physical drive 17 as .a 
function of the remote responses at step 60. As 
10 appreciated by those skilled in the art, the inverse 
model can.be calculated directly from the. input drive 
and remote transducer response. . Also, the term 
"inverse" model as used herein includes a M x N 
"pseudo-inverse " model for . a non-square N x M system. 
5 Techniques for iteratively or directly calculating 
the physical drive 17 to achieve a desired response* of 
the physical system are well known in the art. 

Step 5 6 of applying drive ensembles. to 
physical system 10 can be implemented in accordance 
0 with a first aspect of the present invention discussed 
with reference to FIGS. 5-7. However, this aspect of 
the present invention may be best understood by first 
describing a prior art technique illustrated in FIG . 
4. In FIG. 4; one channel 100 with a collection of 
5 ' drive ensembles is shown. Frames forming the 

collection of drive ensembles, having for example a 
length of five seconds each, are provided sequentially 
as drive 17 input to servo controller 14, and the 
responses 21 are recorded. Illustrated in FIG. 4 are 
0 first, second and third frames for channel 100. 
Typically each frame is made up of a combination of 
periodic waveforms having different frequencies. The 
combinations of periodic waveforms are . repeatable in a 
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subsequent frame . A typical method of measuring the 
FRF of physical system 10 is to generate an orthogonal 
drive where each frame (or ensemble for multiple 
channels) is replicated and the replicated frame is 
5 used for inter- frame smoothing. The FRF analysis is 
then applied to the replicated and unsmoothed frame 
using periodic boundary conditions. Although 
illustrated as one channel, those skilled in the art 
will realize that a frame can represent a collection 
10 of channels in a multi-channel system. 

As described 'above, in the prior art, each 
frame of the channel 100 is repeated. As illustrated 
. in FIG. 4, the. first frame 110 is repeated in frame 
115, the second frame 120 is repeated in frame 125, 
15 and the third frame 130 is repeated in frame 135. To 
avoid discontinuities in the contiguous drive signal,: 
the. first boundary and nearby portions of each new 
non-repeated frame are modified to allow a ramping up 
of the signal and of the response of physical system 
20 10. Thus, as illustrated in FIG. 4, in frames 110, 
120 and 130 respective boundaries 111, 121 and 131 and 
nearby regions of these frames are modified. At 
frames 115, 125 and 135, the first, second and third 
frames are respectively repeated— but in an unmodified 
25 form. Periodic boundary conditions v/hich match frames 
110 and 115 at boundary 112 > frames 120 and 125 at 
boundary 122, and frames 130. and 13 5 at boundary 132, 
allow frames 115, 125 and 135 to remain pure or 
unmodified. Calculating the FRF of the physical 
3 0 system 10 using the relationship between the 
unmodified frames and the response 21 to those 
unmodified frames has been found to be more accurate 
than using modified frames . 
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FIGS. 5-7 illustrate an aspect of the 
present invention in which each frame (ensemble for 
multiple channels) is repeated in an unmodified form- 
at least twice. The additional replications of each 
5 frame allows the system response between identical 
frames to be used to estimate the system repeatable 
bandwidth. In one embodiment illustrated in FIG. 5, 
. each frame of channel 100 is repeated twice. Thus, 
following modified first frame: 110 is unmodified 
10 versions of this frame in frames 115 and 117. 

Consistent with the above description, periodic 
boundary conditions between frames 110 and 115, and 
between frames 115 and 117, allow frames 115 and 117 
to remain unmodified and continuous between adjacent 
15 frames. Similarly, following modified second frame 
120 is unmodified versions of this frame in frames 125 
and 127. Periodic boundary conditions between frames 
120 and 125, and between frames 125 and 127, allow 
frames 125 and 127 to remain unmodified and continuous 
2 0 between adjacent frames. 

FIG. 6 diagrammatically illustrates the 
drive signal of channel 100 in frames 110, 115 and 
117. m As can be seen in FIG. 6, the drive signal in 
frame 110 is modified "near boundary 111 to allow a 
2 5 ramping up of the drive signal and of the system 
response. The periodic boundary conditions allow 
frames 115 and 117 to remain unmodified. 

In a second embodiment illustrated in FIG. 
7 , each frame of the drive signal for channel 100 is 
30 repeated at least three times. Thus, following 
modified first frame ilO is unmodified versions of 
this frame in frames 115, 117 and 118. . Periodic 
boundary conditions between frames" 110 and 115, 
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between frames 115 and 117, arid between frames 117 and 
118, allow frames 115, 117 and. 118 to all remain 
unmodified. One purpose of this . embodiment is to 
illustrate that the. comparison of responses of 
. 5 physical system 10 to different replicants of a 
particular frame or ensemble need not occur between 
adjacent frames. For example, when estimating system 
repeatable bandwidth using the methods described 
below/ the responses to unmodified frames 115 and 118 
10 can be compared. 

The invention is not restricted to 
embodiments in which responses to consecutive 
replicated frames are compared. Although illustrated 
wherein the identical frames comprise a single drive 
15 signal, multiple drive signals each containing one of 
the identical frames can also be used. The multiple 
drive signals can be provided to a single physical 
system at different times, or can be applied to two 
substantially identical physical systems. Reference to 
20 application of identical frames to a physical system 
is intended to cover all of the foregoing. 

Using the methods described above ■ of 
applying drive ensembles to physical system 10, 
system controller 23 estimates the repeatable 

2 5 bandwidth of the physical system. In some 

embodiments, an assumption is made that system non- 
repeatability is only a function of the response 
(i.e., the remote transducers). In these embodiments, 
controller 23 calculates a cross spectral density CSD X 

3 0 from the first unmodified frame or ensemble repeat and 

an additional cross spectral density CSD 2 from the 
second or a subsequent repeat. Cross spectral 

densities (CSDs) are cross-spectrum matrixes between 
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the response signal 21 and the drive signal 17. By 
assuming that the system non-repeatability is only a 
function of the response, the variation AFRF of the 
FRF with respect to' that repeatability may be 
calculated using the relationship shown in Equation 1. 



afrfJ csd >- csd ^ eq 

DSD U 



AFRF represents the lower limit of resolution, of the 
0 FRF and is compared with the FRF which may be 
• calculated as a function of the drive spectral density 
DSD, -which is . the auto- spectrum . matrix of the drive 
signals, using the relationship shown in Equation 2. . 

FRF JCS Dl +CSD 2 ) 

2*DSD Q ' 2 

Any portion of the FRF that lies below AFRF indicates 
that the system has insufficient repeatability. 
Hence, the estimated repeatable bandwidth of the 
0 system should be limited to those regions where the 
FRF is above AFRF. The comparison between FRF and 
AFRF can be implemented 'by controller 23 using the 
following mathematical processes. 

The units scaling from FRF and AFRF are 
removed by factoring out the drive and response units 
using the relationships illustrated in Equations 3 and 
4. 



FRF = R * H * D 



EQ. 3 
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AFRF = R * AH * D EQ. 4 

Matrices R and D are diagonal matrices containing the 
full scales of the drive and response channels, 

5 respectively. Matrices Hand AH are the unitless FRF 
and AFRF matrices. A singular value decomposition of 

//and AH is performed, using the relationships shown 
.in Equations 5 and 6. 

10 H = U * S * . EQ. 5 

AH = 0 * AS * V . EQ. 6 

Where U, U\ V and V are unitary matrices 'and S and 
S are scaling matrices that are diagonal and ordered. 

FIG. 8 is a plot illustrating the first 

15 element 801 of the AS matrix in relation to elements 
805 of the S. matrix for a twelve input . system. FIG. 9 
illustrates the same/ but includes, from the S matrix 
elements 805, ' only the SI and <S12 elements 815 and 
820, respectively. Controller 23 compares the first 

20 element. 801 of the ;AS matrix as a function of 
frequency with the elements 8 05 .of the S matrix. 
Those elements 8 05 of the S matrix that fall below the 
first element 801 of the AS matrix are due to inputs 
that cannot be distinguished by the transducer set 

25 measuring overall system response, and indicate the 
repeatable bandwidth limits of. the system. As 
illustrated at reference number 810, element 820 of 
the S matrix is approximately equal to element 8 01 of 
the AS matrix at a frequency of approximately 38 Hz, 
0 indicating a set of channels of physical system .10 
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that will not be repeatably controllable beyond this 
frequency. If desired, the system repeatable 

bandwidth can be set so that the drive signals 17 on 
all input channels remain below this frequency. In 

5 the alternative, a new model can be constructed by 
eliminating the smallest diagonal (s) of the S matrix 
as a function of frequency, thereby eliminating the 
input frequencies only for the channels of the non- 
repeatable set. 

0 In alternate methods of the present 

invention, an assumption is made that system 
repeatability is a function of only on drive inputs- 17 
to physical system 10. In these embodiments, the 
above-described methods are implemented with the. 

5 variation that AFRF of the FRF is calculated using 
the relationship illustrated in Equation 7. 

AZ7DZ7 (ARSD*CSD'-RSD*ACSD') 

AFRF = y~ : r EQ . 7 

{CSD'*CSD X ) 

0 Where RSD and ARSD are an average and a difference, 
respectively, of the response spectral densities RSD 1 
and RSD 2 of the responses to the first and second 
ensemble repeats. CSD' and ACSD' are an average and a 
difference, respectively, of the -cross-spectral 

5 densities CSD^ and CSD' 2 . 

In accordance with some embodiments of the 
present invention, the FRF measurements are calculated 
with imbedded repeats one channel at a time. Also, in 
some embodiments of the present invention, additional 

) repeats are added to the drive signals and the 
estimate of .the system repeatability is calculated as 
an average of the differences between multiple 
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irepeated drive and response signals. In general, the 
various methods of the present invention can be 
implemented with all channels simultaneously, or one 
channel at a time. 
5 In some embodiments, the environment during 

FRF measurement should be made as similar as possible, 
to the environment of the system during laboratory 
iteration, i.e., in service conditions. This is 
because the largest contributions to system non- 
10 repeatability sometimes come from motions of the 
system itself as compared to an external source. For 
this reason, a differential measurement technique is 
used on a system that is undergoing representative 
multi-channel vibration. Use of this technique can 
15 establish repeatable bandwidth limits of the full 
system. It is possible to alter the model and extend 
the system repeatable bandwidth by eliminating those 
elements of the S matrix that fall below the elements 
of the AS matrix. 
20 In the scope of this repeatability analysis 

it is often desirable to know specifically which 
physical system channels are associated with the 
offending non-repeatable' set of channels as identified 
above. Per the singular value decomposition 

25 H = U * S * V", the system model H for physical 
system 10 is a relationship between drives 17, herein 
denoted as vector X, and responses 21 herein denoted 
as vector Y. Substituting and rearranging 

. H = U * S * V" , Y = H*X 

30 yields 

U'* Y = S * V * X. 
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It has been discovered that matrix V defines which of 
the physical channels X are associated with each 
singular value as function of . frequency. Thus 
specific physical input channels, of a particular set 
5 can be identified, potentially allowing remedial 
action of the input channel (s) to enhance the 
repeatable system bandwidth. For instance/ with 
reference to FIG. 1, remedial action could include 
hardware fixes of the actuator (s) 15 or interfaces 16, 
10 corrections of problems with the servo controller 14, 
correction of electrical problems between the system 
controller 23 and the .servo controller 14,. to name a 
few. . 

When all input channels are determined to be 
15 working properly, the problem can be due to the 
response transducer configuration. In this case, the 
column of; V associated with a small singular value 
defines a combination of input channels, possibly at' a. 
high amplitude, that is not detected by the response 
20 transducer set. In other words, all of the response 
transducers miss this state of vibration. 

Knowing the input combination that produces this 
effect, can lead directly to a corrective modification 
of the response transducer set. In difficult cases, 
25 the physical system can be driven with such a 
combination of inputs to: produce an observable level 
of vibration> but insignificant response '. from the 
transducers. The location of one or more significant 
response points might be determined by observation, 
30 greatly accelerating the process of adding or moving 
transducers to achieve an acceptable set . 

Although, steps 56 and 57 can be performed 
during system identification 55, as described above, 
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where the drive ensembles are designed to provide 
information required to calculate the system model. . 
However/ the drive ensembles can be those used during 
iterative phase 60 where the drive ensembles are 
5 calculated during the iterative process to reproduce 
the desired response. The . drive ensembles from the 
iterative phase 60 would not include repeated 
sequences of modified and unmodified frames, but 
rather are the entire drive signal. As those skilled 
in the art would appreciate, the iteration, drive 
ensembles may not include broad band frequency 
information that is inherent to the drive ensembles 
that are used for system modeling . . 

Although . described above wherein an 
estimation of system repeatability to the differences 
in the responses to substantially identical drive 
ensembles was described using the singular value 
techniques presented, other techniques for inferring" 
or estimating repeatable system bandwidth can be 
utilized in accordance with the invention. For 
example, these techniques include, but are not limited 
to, frequency domain analysis such as AFRF, ACSD, 
ARSD, and AS. of the responses to the identical drive 
ensembles; statistical domain analysis such as ARMS 
(root mean square of the response time histories) of 
the responses to the identical drive ensembles; and 
time history domain analysis such as the difference 
between the responses to the identical drive 
ensembles , 

Although the present invention has been 
described with reference . to preferred embodiments, , 
workers skilled in the art will recognize that, changes 
may be made in form and detail without departing from 
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the spirit and scope of the invention. For example, 
those of skill in the art will recognize that a slight 
modification to the "unmodif ied" drive ensembles will 
.still result in the benefits of* the invention. 
5 Further, the application of identical" drive 

ensembles for use in calculating or estimating a 
repeatable system bandwidth can implemented with 
substantially or approximately identical drive 
ensembles in alternate embodiments. Further still, 
10 application of the identical drive ensembles to the 
•physical system at different times can be implemented 
using any technique (e.g. in separate drive signals or 
in appended drive signals) . . Also, the various methods 
of the present invention, can be implemented with all 
15 channels simultaneously, or one channel at a time. 
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WHAT IS CLAIMED IS: 

1. A method of identifying characteristics of a 
physical system, the method comprising: 

applying substantially identical drive ensembles to 
the physical system and obtaining corresponding 
responses from the physical system; and 

estimating a repeatable bandwidth of the physical 
system as a function of the applied drive 
ensembles and the corresponding obtained 
responses. 

2. The method of claim 1/ and further comprising 
calculating simulation drive inputs for the physical 
system . as a function of the estimated repeatable 
bandwidth. of the physical system in order to obtain a 
desired response from the physical system. 

3. The method of claim 1, wherein applying the 
substantially identical drive ensembles to the 
physical system further comprises applying a modified 
drive ensemble to the physical system prior to 
applying the first and, second substantially identical 
drive ensembles to the physical system, wherein the 
first and second substantially identical drive 
ensembles are substantially unmodified versions of the 
modified drive ensemble. 

4. The method of claim 3, wherein applying the first 
and second substantially identical drive ensembles 
further comprises applying the second drive ensemble 
immediately after applying the first drive ensemble. 
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5. The method of claim 3, wherein applying the first 
and second substantially identical drive ensembles 
further • comprises applying a third drive ensemble 
between the first and second substantially identical 
drive ensembles. 

6. The method of claim 3, wherein applying the. first 
and second substantially identical drive ensembles 
further comprises : 

applying a first drive signal having the first drive 
ensemble; and 

a PPlyi n 9" a second drive ensemble having the second 
drive ensemble . 

7. The method of claim 1, wherein estimating a 
repeatable bandwidth includes at least * one of 
frequency domain analysis, statistical domain analysis 
and time history domain analysis of the responses to 
the substantially identical drive ensembles. 

8. The method of claim 7, wherein estimating the 
repeatable bandwidth of the physical system further 
comprises 

calculating a first ' cross spectral density CSD X 

between the' first drive ensemble and the 

corresponding response; 
calculating a second cross spectral density CSD 2 

between the second drive ensemble and the 

corresponding response; and 

calculating a variation AFRF of a frequency response 
function FRF for the physical system as a 
function of the first cross spectral density S1SD 1 
and the second cross spectral density CSD 2 . 
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9, The method of claim 8, wherein the frequency 
response function FRF is calculated as a function of 
the first cross spectral - density CSD lt the second 
cross spectral density CSD 2 , and of a drive spectral 
density DSD, using the relationship: 

FRF _ (CSD l +CSD 2 ) 
2 * DSD 

10 v The method of claim 8, wherein calculating the 
variation AFRF of the frequency response function FRF 
for the physical system further comprises calculating 
the variation AFRF as a function of the first cross 
spectral density CSD X , the second cross spectral 
density CSD 2 , and a drive spectral density DSD using 
the relationship: 

(CSD l ^CSD 2 ) 



AFRF = ■ 



DSD 



11. The method of claim 8/ wherein calculating the 
variation AFRF of the frequency response function FRF 
for the physical system further comprises calculating 
the variation AFRF as a function of the first cross 
spectral density CSD X , the second cross spectral 
density CSD 2 , and a drive spectral density DSD using 
the relationship: 

._■„ (ARSD*CSD'-RSD*ACSD') 
At Kr — 7 c 

(CSD'*CSD') 

where RSD is an average of the response spectral 
densities RSD! and RSD 2 of the responses to the first 
and second substantially identical drive ensembles, 

wherein ARSD. is a difference of the response spectral 
densities RSD X ancl RSD 2 of the responses to the first 
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and second substantially identical drive ensembles, 
CSD' is an average of the cross-spectral densities 
CSDV and CSD' 2 of the responses to the first and 
second substantially identical drive ensembles, and 
ACSD' is a difference of the cross-spectral densities 
CSD', and CSD' 2 of the responses to the first and 
second substantially identical drive ensembles. 

12. The method of claim 8, wherein estimating the 
repeatable bandwidth for the physical system further 
comprises comparing a singular value decomposition as 
a function of frequency of the variation AFRF to a 
singular value decomposition as a function of 
frequency of the frequency response function FRF. 

13. The method of claim 12 , and further comprising 
identifying inputs to the physical system to which the 
physical system is insensitive as a function of 
frequency using the relationship: 

U' * Y = S * V * X 
where H is a model of the physical system, X is a 
matrix, of inputs of the physical system, Y is a matrix 
of responses of the physical system such that Y=H*X 

and H = U * S * V", . and where U' and V are 
rotational matrices, V' being indicative of the 
inputs to the physical system to which the physical 
system is insensitive. . 

14. A computer readable medium, including instructions 
readable by a computer, which when implemented, cause 
the computer to identify characteristics of a physical 
system, the instructions performing steps comprising: 
applying substantially identical drive ensembles to 
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the physical system and obtaining corresponding 
responses from the physical system; and 
estimating a repeatable bandwidth of the physical 
system as a function., of the applied drive 
ensembles and • the corresponding obtained 
responses. 

15. The computer readable medium of claim 14, and 
further including instructions performing step 
comprising calculating simulation drive inputs for the 
physical . system as a function of the estimated 
repeatable bandwidth of the physical system in order, 
to obtain a. desired response from the physical system. 

16. The computer readable medium of claim 14, wherein 
the instructions performing step of applying the 
substantially identical drive* ensembles to the 
physical system further comprises applying a modified 
drive ensemble to the physical system prior to 
applying the first and second substantially identical 

•drive ensembles to the physical system, wherein the 
. first and second substantially identical drive 
ensembles are substantially unmodified versions of the 
modified drive ensemble. y ' 

17. The computer readable medium of claim 16, wherein 
the instructions performing step of applying the first 
and second substantially identical drive ensembles 
further comprises applying the second drive ensemble 
immediately after applying the first drive ensemble . 

18. The computer readable medium of claim 16, wherein 
the instructions performing step . of applying the first 
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and second substantially identical drive ensembles 
further comprises applying a third drive ensemble 
between the first and second substantially identical 
drive ensembles. 

19.. The computer readable medium of claim 16, wherein 
the instructions performing step of applying the first 
and second substantially identical drive . ensembles 
further comprises: 

applying a first drive signal having the first drive 
ensemble; and 

applying" a second drive ensemble having the second 
drive ensemble. ' 

20. The computer readable medium of claim 14, wherein 
the instructions performing "step of estimating a 
repeatable bandwidth includes at least one ■ of 
frequency domain analysis, statistical domain analysis 
and time history domain analysis of the responses to 
the substantially identical drive ensembles. 

21. The computer readable medium of claim 20, wherein 
the instructions performing step of estimating the 
repeatable bandwidth " of the physical system further 
comprises : 

calculating a first cross spectral density ■CSD 1 

between the. first drive- ensemble and the 

corresponding response; 
calculating a second cross spectral density CSD 2 

between the second drive ensemble and the 

corresponding response; and 

calculating a variation AFRF of a frequency response 
function. FRF for the physical system as a 
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function of the first cross spectral density CSD 1 
and the second cross spectral density CSD 2 . 

22. The computer readable medium of * claim 21, wherein 
the frequency response function FRF is calculated as a 
function of the first cross spectral density CSD 1 , the 
second cross spectral density CSD 2 , and of a drive 
spectral density DSD, using the relationship: 

FRF (CSD, + CSD 2 ) 
2* DSD 

23. The computer readable medium of claim 21, wherein 
the step of calculating the variation AFRF of the 
frequency response function FRF for the physical' 
system further comprises calculating the variation 
AFRF as a function of the first cross spectral density 

' CSD 1§ the second cross spectral density CSD 2 , and a 
drive spectral density DSD using the relationship: 

^ RF JCSD X -CSD 2 ) 
DSD 

24. The computer readable medium of claim 21,. wherein 
the step of calculating the variation AFRF of the 
frequency response function FRF for .the physical 
system further comprises calculating the variation 
AFRF as a function of the first cross spectral density 
CSD A , the second cross spectral density CSD 2 , and a 
drive spectral density DSD using the relationship: 

£j? RF = (^SD * CSD'-RSD* ACSD') 
(CSD'*CSD') 

where RSD is an average of the response spectral 
densities RSD 2 and RSD 2 of the responses to the first 
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and second substantially identical drive ensembles, 
wherein ARSD is a difference of the response spectral 
densities RSD X and RSD 2 of the responses to the first 
and second substantially identical drive enssmbles, 
.CSD' is an average of the cross-spectral densities 
CSD' X and CSD' 2 of the responses to the first and 
second substantially identical drive ensembles, and 
ACSD' is a difference of the cross-spectral densities 
CSD', and CSD' 2 ' of the responses to the first and 
second substantially identical drive ensembles. 

25. - The computer readable medium of claim 21, wherein 
the step of estimating the repeatable bandwidth for 
the physical system further comprises comparing a 
singular value decomposition as a function of 
frequency of. the variation AFRF to a singular value 
decomposition as a function of frequency of the 
frequency response function FRF. 

26. the computer readable medium of claim 25, and 
further comprising an instructions performing step * of 
identifying inputs to the physical system to which the 
physical system is insensitive as a function of 
frequency using the relationship: 

U' * Y = S * .V * X 
where H is a model of the physical system, X is a 
matrix of inputs of the physical system, Y is a matrix 
of responses of the physical system such that Y=H*X 

and H = U * S * V", and where U' and V are 
rotational matrices, V being indicative of the 
inputs to the physical system to which the physical 
system is insensitive. 
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27. An identification system for identifying 
characteristics of a physical system, the 
identification system comprising: 

an input coupled to the physical system and adapted to 
receive responses from the physical system ,- 

an output coupled to the physical system and adapted 
to provide drive ensembles to the physical 
system; and 

a controller coupled ■ to the input and to the output 
and adapted to apply substantially identical 
drive ensembles to the physical system, to obtain 
corresponding responses from the physical system, 
- and to estimate a repeatable bandwidth of the 
physical system as a function of the applied 
drive ensembles and the corresponding obtained 
responses . . . 

28. An identification ' system for identifying 
characteristics of a physical system, the 
identification /system comprising: 

input means for receiving responses from the physiical 
system; 

output means for providing drive ensembles to the 
physical system; and 

control means for applying substantially identical 
drive ensembles to the physical system, obtaining 
corresponding responses from the physical system, 
and estimating a repeatable bandwidth of the 
physical system as a function of the applied 
drive ensembles and the corresponding obtained 
responses . 
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